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I ABSTRACT 

, Investigation of the behaviour of accreting matter close to the black hole event horizon 

is of fundamental importance in relativistic and high energy astrophysics because they 
^ I provide the key features of the diagnostic spectra of the stellar mass and super-massive 

^ . black holes. In this paper, we examine the terminal behaviour of general relativis- 

' tic matter in multi-transonic, advective black hole accretion discs. We compute, for 

I the first time we believe, the values of various dynamical and thermodynamic multi- 

1/^ . transonic flow variables extremely close (< 0.01 rg) to the event horizon and study the 

CNJ ' dependence of these variables on fundamental accretion parameters. Our calculation 

^ is useful for a better understanding of Hawking radiation from acoustic black holes. 

m . 

. Key words: Accretion, accretion discs — black hole physics — general relativity — 

' hydrodynamics — Hawking Radiation 
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C/3 1 INTRODUCTION 

IL* Gravitational capture of surrounding fluid by massive astrophysical objects is known as accretion. There remains a major dif- 
. 5^ ' ference between black hole (BH) accretion and accretion onto other cosmic objects including neutron stars and white dwarfs. 
, For celestial bodies other than black holes, infall of matter terminates either by a direct collision with the hard surface of the 

5_j accretor or with the outer boundary of the magneto- sphere, resulting the luminosity through energy release from the surface. 
Whereas for black hole accretion, matter ultimately dives through the event horizon from where radiation is prohibited to 
escape according to the rule of classical general relativity (GR) and the emergence of luminosity occurs on the way towards 
the black hole event horizon. The efficiency of accretion process may be thought as a measure of the fractional conversion of 
gravitational binding energy of matter to the emergent radiation and is considerably high for black hole accretion compared 
to accretion onto any other astrophysical objects. Hence accretion onto classical astrophysical black holes has been recognized 
as a fundamental phenomena of increasing importance in relativistic and high energy astrophysics (Frank, King & Raine 
1992, hereafter FKR, Shapiro & Teukolsky 1983). The extraction of gravitational energy from the black hole accretion is 
believed to power the energy generation mechanism of X-ray binaries and of the most luminous objects of the Universe, the 
Quasars and active galactic nuclei (AGN). The BH accretion is, thus, the most appealing way through which the all pervading 
power of gravity is explicitly manifested. If the infalling matter does not possess intrinsic angular momentum, accretion flow 
remains spherically symmetric. However, Inter stellar / intergalactic fluid is always likely to posses non-vanishing rotational 
energy, sufficient to dynamically break the spherical symmetry, and in almost all real physical situations, accreting matter is 
thrown into circular orbits around the central accretor, leading to the formation of the accretion disc around the galactic and 
extra-galactic black holes. 

If the instantaneous dynamical velocity and local acoustic velocity of the accreting fluid, moving along a space curve parame- 
terized by r, are u{r) and a{r) respectively, then the local Mach number M{r) of the fluid can be defined as M(r) = The 
flow will be locally subsonic or supersonic according to M{r) < 1 or > 1, i.e., according to u(r) < a(r) or it(r) > a(r). The 
flow is transonic if at any moment it crosses M — 1. This happens when a subsonic to supersonic or supersonic to subsonic 
transition takes place either continuously or discontinuously. The point (s) where such crossing takes place continuously is 
(are) called sonic point(s), and where such crossing takes place discontinuously are called shocks or discontinuities. In order 
to satisfy the inner boundary conditions imposed by the event horizon, accretion onto black holes exhibit transonic properties 
in general, which further indicates that formation of shock waves are possible in astrophysical fluid flows onto galactic and 
extra-galactic black holes (Das, Pendharkar & Mitra 2003, and references therein. Das 2001, and references therein). The study 
of dynamical behaviour of such transonic accretion flow near the BH event horizon is considered to be immensely important. 
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Due to the strong curvature of space-time close to the black hole, accreting fluid is expected to show extreme behaviour just 
before plunging into the event horizon; this tremendously hot, ultra-fast flow with it's high-density proflle is supposed to 
provide the key features of the diagnostic high energy spectra of galactic and extra-galactic black hole candidates. 
Following the pioneering contributions of Bardeen, Press & Teukolsky (1972) and Novikov & Thorne (1973, NT hereafter), a 
number of papers deal with the general rclativistic (GR) BH accretion disc. The initial attempt by Fukuc (1987 and references 
therein) and Chakrabarti (1990, 1996, hereafter C96) to analytically study the GR transonic BH accretion was followed by 
other independent analytical works (Kafatos & Yang 1994, Yang & Kafatos 1995, YK hereafter, Pariev 1996, Peitz & Appl 
1997, Lasota & Abramowicz 1997, hereafter LA, Lu, Yu, Yuan & Young 1997). Although they made a number of important 
contributions to the study of general relativistic advective black hole accretion discs, none of the above mentioned works 
investigates the flow behaviour close to the event horizon in detail. Gammie and Popham (Gammie & Popham 1998, Popham 
& Gammie 1998) performed the only self-consistent analytical work present in the literature so far which exploits the complete 
general relativistic formalism to study various flow variables near the event horizon. Recent interesting work of Becker & Le 
(2003) uses post-Newtonian asymptotic analysis to study the properties of inner region of advective accretion flows, and their 
results are in good agreement with general relativistic description. However, neither Gammie ot al nor Becker ct. al. captures 
the multi-transonic properties of the flow which is of considerable astrophysical importance. Only the multi-transonic flow can 
produce shock waves in post-Newtonian black hole accretion discs (Das, Pendharkar & Mitra 2003, and references therein). 
Shock waves in rotating flows may provide an efficient mechanism for conversion of significant amount of the gravitational 
energy into radiation by randomizing the directed infall motion of the accreting fluid and the hot and dense post-shock flow is 
considered to be a powerful tool in understanding various important astrophysical phenomena like the spectral properties of 
the black hole candidates (Chakrabarti & Titarchuk 1995, and references therein), the formation and dynamics of accretion 
powered cosmic jets (Das & Chakrabaxti 1999, and references therein. Das, Rao & Vadawale 2003), the origin of QPOs in 
galactic sources (Das 2003, and references therein). One thus understands the pressing needs for a self-consistent analytical 
model capable of studying the behaviour of general relativistic, multi-transonic accretion flow sufficiently close to the event 
horizon. 

Motivated by the above mentioned arguments, in this paper we concentrate on the stationary, axi-symmetric, complete gen- 
eral relativistic accretion solutions in Schwarzschild metric which contains multiple critical points, and for such solutions, 
we calculate all relevant dynamical and thermodynamic flow variables extremely close to the event horizon. Such a solution 
scheme for purely spherical general relativistic mono-transonic BH accretion was outlined in Das 2002. We also exhaustively 
study the dependence of such variables (close to the event horizon) on all important initial boundary conditions governing 
the flow. Our generalized calculation is valid for both super- as well as for sub-Eddington accretion onto BH of any mass. 
We thus, provide a useful and self-consistent procedure to study the terminal behaviour of matter in multi-transonic, general 
relativistic advective accretion discs around BHs, which has never been done in any of the existing works on black hole accre- 
tion discs. Further calculations, with an ultimate ambition towards modelling the most general viscous, transonic, shocked 
hydromagnetic accretion in Kerr-Neumann space time is in progress (T. K. Das, P. J. Wiita & P. Barai, in preparation) and 
would be reported elsewhere. 

At this point, we would like to mention that instead of dealing with the so called standard disc model with signiflcantly 
high angular momentum (where the presence of viscous stress allows the infall of accreting material onto the black hole, i.e., 
outward viscous transport of angular momentum takes place to weaken the centrifugal barrier), we rather mainly concentrate 
on accretion flows with relatively low intrinsic angular momentum (sub-Keplerian angular momentum distribution) where 
substantially significant advection velocities may be obtained even for practically inviscid flow, although our formalism allows 
sufficient flexibility to incorporate the viscous transonic flows as well, see §4.1 for more detail. Such weakly rotating flows have 
not been explored much in the literature, although they are well exhibited in nature for various real physical situations like de- 
tached binary systems fed by accretion from OB stellar winds (lUarionov & Sunyaev 1975; Liang & Nolan 1984), semi-detached 
low-mass non-magnetic binaries (Bisikalo et al. 1998) and supermassive BHs fed by accretion from slowly rotating central 
stellar clusters (Illarionov 1988; Ho 1999 and references therein). Even for a standerd Keplarian accretion disc, turbulence 
may produce such low angular momentum flow (see, e.g., Igumenshchev & Abramowicz 1999, and references therein). 



2 FORMALISM 

2.1 The energy momentum tensor 

To provide the most general description of fiuid fiow in strong gravity, one needs to solve the equations of motion for the 

fluid and the Einstein equations. The problem may be made simplified by assuming the accretion to be non-self gravitating 
so that the fluid dynamics may be dealt in a metric without back-reactions. Unless otherwise mentioned for any specific 
reason(s), hereafter we define Schwarzschild radius rg = ^S^^SiL^ where Mbh is the mass of the black hole, G is the Universal 
gravitational constant and c is the velocity of light. The radial distances and velocities are scaled in units of rg and c respectively 
and all other derived quantities are scaled accordingly; G = c = Mbh = 1 is used. We use the Boyer-Lindquist co-ordinate 

with signature h H — h, and an azimuthally Lorentz boosted orthonormal tetrad basis co-rotating with the accreting fluid. 

We define A to be the specific angular momentum of the flow and neglect any gravo-magneto-viscous non-alignment between 
A and BH spin angular momentum. 

Let be the four velocity of the (perfect) accreting fluid. The energy momentum tensor 9'"' of such a flow could be written 
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Q'^" = (p + p) Vf^v,, + pg^v, or T = (p + p) v®v + pg (la) 

A complete description ol flow behaviour could be obtained by taking the co- variant derivative of Ss*^" and pv'^ to obtain the 
energy momentum conservation equations and the conservation of baryonic mass. 

However, at this stage, the complete solution remains analytically untenable unless we axe forced to adopt a number of simpli- 
fied approximations. In this paper, as already mentioned in §1, we would like to study the inviscid accretion of hydrodynamic 
fluid in Schwarzchild metric. Hence, our calculation will be focused on the stationary axisymmetric solution of the following 
equations: 



= 0, (pt;")^ 



= (16) 



where a is the Kerr parameter, and the semicolon (;) denotes the co-variant derivative. The radial momentum balance 
condition may be obtained by making the expression [v^v,, + g^v) ^^t" to be equal to zero. The 'stationarity' condition 
imphcs the vanishing of the temporal derivative of any scalar field in the disc or the vanishing of the Lie Derivative of any 
vector or tensor field along the killing vector e.g., ^ = if 11 is a scalar field, and Zla/atll = if n represents a vector 
or a tensor field. The 'axisymmetry' property is endowed with the space like Killing field {-§^Y ■ 

Exact Solution of the above mentioned conservation equations requires the knowledge of accretion geometry, as well as 
introduction of a suitable equation of state. We concentrate on polytropic accretion for which p = Kp'' , where K, p and 7 are 
the rnonotonic and continuous function of the specific entropy density, the adiabatic index and the rest mass density of the 
fiow, respectively. However, one may note that the polytropic accretion is not the only choice to describe the general relativistic 
transonic black hole accretion and equations of state other than the adiabatic one has also been used by several authors, like 
YK (isothermal equation) or Manmoto 2000 (two temperature plasma). There is no reason that more realistic equation of 
state can not be employed except the fact that the analysis would be algebraically more complicated and hence would be 
dominated by additional numerical details. Specific proper flow enthalpy is taken to be /i = (7 — 1) {7 — (l -|- o^) } , where 
a is the polytropic sound speed defined as (Weinberg 1972, NT, Frank, King & Raine 1992) a = [dp/de^J'^ = ^'i (r(r),7) = 
*2 ip,p,7), where T(r) is the local fiow temperature, e is the mass-energy density, and ^"2} are known functions. The 
subscript <S indicates that the derivative is taken at constant specific entropy. 



2.2 The disc height 

We assume that the disc has a radius dependent local thickness H(r), and it's central plane coincides with the equatorial plane 
of the BH. It is a standard practice (Matsumoto et. al. 1984, Paczy'nski 1987, Abramowicz, Czerny, Lasota & Szuszkiewicz 
1988, Chen & Taam 1993, KY, Artemova, Bjornsson & Novikov 1996, Wiita 1999, Hawley & Krolik 2001, Armitage, Reynolds 
& Chiang 2001) to use the vertically integrated model in describing the black hole accretion discs where the equations of 
motion apply to the equatorial plane of the BH, assuming the flow to be in hydrostatic equilibrium in transverse direction. 
We follow the same procedure here. Flow variables, if a vector or tensor, are Lie dragged along (^) while averaging over 
all 4>, and then are averaged over a length scale of the order of H(r). We consider the flow to be 'advective', i.e., to possess 
considerable radial three-velocity. 

Following the process outlined in LA, we calculate the disc height. In Newtonian framework, the disc height in vertical 
equilibrium is obtained from the z component of the non-relativistic Euler equation where all the terms involving velocities 
are neglected and the higher powers of are ignored in the Taylor's expansion. For general relativistic disc, the vertical 
pressure gradient in the co-moving frame is compensated by the tidal gravitational field. We assign the Lorentz factor F to 
be equal to: 



in the Schwarzschild metric and approximate || — » |, so that the z component of the general relativistic Euler equation 
provides: 



Hir) = 2V2^ P^'-'-^^'- + ^"^ (26) 

It is trivial to show that for the equation of state and the metric used in this paper, 

_ pa^ (7 - 1) . . 

P-y_^(l + a2) ^^^^ 

Using the above equation and the relation (NT, FKR) o^(r) = ''^^^^ = 0^T(r), where where = ^J'-^^^ M is the mean 
molecular weight, ruH'-^mp is the mass of the hydrogen atom and k is Boltzmann's constant, we finally obtain the disc height 
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H{r) = 5.66eT 2 (r) 



^){ 



(r-1) 



7- (l + e2T(r)) 



(2d) 



2.3 The conserved specific energy of accretion flow 



Temporal (zeroth) component of the first part of eq. (lb) leads to the conservation of specific flow energy £ (relativistic 
analogue of Bernoulli's constant) along each streamline as (Anderson 1989) £ = hvt, where is the four velocity. Prom the 
normalization condition of the four velocity {VfiV^ = — 1), one writes (NT): 

2 ^ 4r (r - 1) 

* (1 - w2) (1 - 2QX) {g^^ + 2\gt4,) ^ ' 

where H (= «**/«') is the angular velocity and u is the radial three velocity in the co-rotating fluid frame. 
One calculates the required metric elements as: 

r — \ 2 

gtt = , g4,<i, = 4r , gt4, = g,i,t = (4a) 

r 

to obtain the expression for the angular velocity as: 

Hence vt for a non-spinning (Schwarzschild) black hole may be derived as: 



Vt 



Schwarzschild 



2GMbh 
7i 



''^/(l-w2)(r3 + A2- AV) ^^''^ 



Using the value of specific fiow enthalpy h (defined in §2.1), eq. (4c), and u = aM(r) = QT^ (r)M(r), M(r) being the radial 
flow Mach number, we obtain the expression for the conserved specific energy (which includes the rest mass energy) as: 



c_ ^(7-1) / r_-l 

[7 - {1 -h e2r(r)}] V {1 - Q^T{r)M^{r)} {r^ - AV + A^) 



2.4 The mass and entropy accretion rate 

We obtain the mass accretion rate Min by integrating the second part of the eq. (lb): 



- 1) ■ 

1-m2(j.)' 



Min = 47r«(r)^/:p^-3^S(r)^(^) (5a) 



where S(r)^(^j is the vertically integrated surface density of the disc material. Hence one writes: 

nr)Hir) = / Pir)dz (56) 

J -H(r)/2 

where p(r) in is the mean value of radial density on the equatorial (z = 0) plane. We use the value of H{r) from eq. (2c) to 
perform the above integration and then substitute the value of S(r)^(^^ in eq. (5a), and use u{r) = a{r)M{r) = QT'^M{r) 
to flnally obtain the mass accretion rate as: 



The entropy accretion rate H can be defined as a quasi-constant multiple of the mass accretion rate: 

■k= p^p^ M,n =KT^ Min (5d) 

After substituting the value of Min from eq. (5c) with the transformation T{r) — > "q'^ , M{r) = {w(r),a(r)}, H reads: 
^ „^ 7+1 /r (r - 1) {7-3 - A2 (r - 1)} f 0^(7-1) ] 

S = 71.057^«^^ ^ | ,-(] + j) | (5^) 

Note that, in absence of creation and annihilation of matter, while the mass accretion rate is an absolute constant of motion, 
entropy accretion rate is not. As the expression for H contains K (= pp''), which is a measure of the entropy density of the 
fiow, H remains constant throughout the flow only if the local entropy density of the flow remains unchanged. Thus, H is a 
constant of motion for shock-firee polytropic accretion and becomes discontinuous (increases) at the shock location, if shock 
forms in the accretion. Thus for a shock free non-dissipative flow, ^ = along a streamline. 
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2.5 Dynamical velocity gradient and sonic quantities 

We take the logarithmic differentiation (with respect to r) of both the side of eq. (5e) to obtain: 



da{r) 
dr 



7 



du 



where: 



u (u^ — 1) dr 



A2 - 3r2 



+ f{r,X) 



2r (r - 1) 2 (r^ - Xh- + X^) 



(6a) 



(66) 



We now transform T{r) — > air), M{r) — > {u{r),a{r)} in the expression for 8, differentiate both side of eq. (4d) and 
substitute the value of da{r) /dr from eq. (6a) to obtain the dynamical velocity gradient as: 



2r 



ril + A^Cl-r) 



2r-l _ 2a^ 
2r(r-l) 7 + 1 



2r(r-l) + 2[r3 + A2{l-r)] 



V 



u(u2-l)(7+l) 



+ 



(6c) 



Hereafter, we use the notation [Pa] for a set of values of {f , A, 7}, which will be our three parameter initial boundary condition 
determining the flow behaviour. Since the flow is assumed to be smooth everywhere, if the denominator of eq. (6c) vanishes at 

any radial distance r, the numerator must also vanish there to maintain the continuity of the flow. One therefore arrives at the 
so called 'sonic point (alternately, the 'critical point') conditions' by simultaneously making the numerator and denominator 
of eq. (6c) equal zero. The sonic point conditions can be expressed as: 



7+1 



1 


2r3-A2 


2r - 1 


2r 


ra + A^{l — 1) 


2r{r-l) 



(r-l) + 2[r3+A2(l-r)J 



(6d) 



where the subscript s indicates that the quantities are to be measured at the sonic point (s). For a fixed [Ps], we substitute 
the values of as and in the expression of £ with the substitution T(r) = o^(r)/0^, M(r) = {^(r), a(r)}, and obtain the 
following polynomial, solution of which provides the sonic point (s) rs'- 



where 



and 



r?(rs,A) = 



1 + 



(rs,A) 
7+1 



Ts (7 - 1) 



7-r7(rs,A) 
*(rs,A) 







(6e) 



*(rs,A) = 



1 






2r,, -1 


2rs 


|_r|+A2(l-r)J 




2r,, (r,, -1) 



l-2r3 , A2-3rg 
2r-(ra-l) 2[r|+A2(l-r-3)J 



To determine the behaviour of the solution near the sonic point, one needs to evaluate the value of (^) at that point (the 
'critical velocity gradient' (^) ) by applying L 'Hospitals' rule to eq. (6c): 



rdu\ 



m 



(6/) 



We use the expression of A/" and V from eq. (6c), calculate Ltr^r, {dAf /dr) and Ltr^rs {dT>/dr). After a number of complicated 
algebraic manipulations, eq. (6f) reduces to the following polynomial which could be solved to obtain the critical velocity 
gradient: 



2(27-30?) /d„N2 

~2[^] +4C(rs,A) 



(7+l)(Mi-l)' \dr 



Mi - 1 



rdu\ 



+ 



^ Os^(rs,A) 



7+1 



2?(rs,A) 



1 + a. 



7+1 



2rs 



2.ri - X' 



rs {rs - 1) rf + A^ (1 - rs) 



2Qri 



12rl - 2)^ (2,rs - 2) 



5rt 



4rf - A2 (3r| - ^rs + 1) 



= 



(6ff) 
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where 



?(rs,A) = 



1 - 2r - 3r^ 



2r (r - 1) 2 [r^ + A2 (1 - r)] 



3 RESULTS 

3.1 Parameter space for multi-transonic accretion 

For a particular [Pa], if ^ [Pa] denotes the universal sot representing the entire parameter space covering all possible val- 
ues of [Pa], and if S [Pa] represents one particular subset of .4 [Pa] which contains only the values of [Pa] providing more 
than one real roots of eq. (6e), then BlVs] can further be decomposed into two subsets C ['Pa] and 'DlVs] such that 



ClVs] C B[V-^] for H(r,„) > S{rout), and ©[Pa] Q B[Vz] for H (r,„) < H (ro„t) , then for [Pa] G C [Pa], we get 



multi-trans omc accretion where two real physical inner and outer (with respect to the BH location) X type sonic points Vin 
and Tout encompass one O type unphysical middle sonic point Vmid in between. One can obtain multi-transonic winds for 
S (rin) < H (rout), e.g., [Pa] G V [Pa] provides multi-transonic winds. 

In Figure 1, the wedge-shaped surfaces represent (f ,A) G [Pa] G C [Pa] C B [Pa] for three different 7 = 1.33 (area bounded 

by dotted lines), 1.43 (area bounded by dashed lines), 1.53 (area bounded by solid lines). If Smax be the maximum value of 
the energy and if Xmax and Xmin be the maximum and minimum values of the angular momentum respectively for C [Pa] for 
a fixed value of 7, then [£max,Xmax,Xmin] non-linearly anti-correlates with 7. It is obvious from the figure that, as the fiow 
makes a transition from its ultra-relativistic * to its purely non-relativistic limit, the area representing C [Pa] decreases. Using 
the formalism developed in §2, it is easy to obtain the similar regions in [Pa] space representing the multi-transonic winds 
as well. However, as our main interest is to study the behaviour of the mfalling matter close to the black hole, we do not 
concentrate on wind solutions in this paper and hence do not provide the representative surfaces for [Pa] GV [Pa]. 

3.2 Integral curves of motion 

In Figure 2, we show the integral curves of motion for multi-transonic accretion. While the distance from the event horizon of 
the central BH (scaled in the units of Vg and plotted in logarithmic scale) is plotted along the X axis, the local Mach number 
of the fiow is plotted along the Y axis. For values of C [Pa] shown in the figure, ABCD represents the accretion passing 
through the outer soruc point B, location of which can be found by solving eq. (6e). EBI represents the self-wind. Flow along 
GFH passes through the inner sonic point F and encompasses a middle sonic point rmid location of which is shown in the 
figure using an asterisk. Similar topologies could be obtained for any other [Pa] G C [Pa]. The overall scheme for obtaining the 
above mentioned integral curves is as follows: 

First we compute rin, rmid and Vout by solving eq. (6e). Then we obtain the dynamical velocity gradient of the flow at sonic 
points by solving eq. (6g). For a chosen A4in (scaled in the units of the Eddington rate MEdd), we then compute the local 
dynamical fiow velocity tt(r), the local polytropic sound speed a(r), the local radial Mach number M(r), the local fiuid density 
p(r) and any other related dynamical or thermodynamic quantities by solving the eq. (6a-6g) from the outer as well as from the 
inner sonic point using fourth order Runge-Kutta method. Flows passing through rout may generate AH = H (rin) — 3 ij'out) 
amount of entropy density by undergoing a steady, standing, Rankiiie-Hugoniot type of shock transition to produce subsonic, 
hotter and shock compressed post-shock solutions (branch GFH) which dives onto the event horizon supersonically after 
passing through nn- Such shocks are time- like three-surfaces of first-order discontinuities, and axe formed if the following 
equation is satisfied: 



where E(= M_/M+) and H(= S+/S_) are shock compression and the entropy enhancement ratio at the shock and T{— /+) 

and u(— /+) arc the prc-/post-shock temperature and dynamical velocities of the fiow respectively. If [Pa] G S [Pa] provides 
real solutions of eq. (7), it directly foUows that S [Pa] C C [Pa] with {[Pa] G N [Pa] | [Pa] G C [Pa] , [Pa] i S [Pa]} representing 
the region producing non-stationary shock solutions. 

3.3 Computation of the quasi-terminal values 

This section refiects the main findings of our paper. We define 'quasi-terminal values' (QTV) to be the value of any fiow 
variable Vf as V^''"^ at a distance rs = rg (1 + S) with < 5 « 1. Terminal value {S = 0) of Vf diverges due to singularity 
at event horizon. For any [Pa] G C [Pa], we provide simultaneous numerical solution of eq. (4d - 6g) to calculate various V^^^ . 
The basic scheme of such calculation is the following: 

* By the term 'ultra-relativistic' and 'purely non-relativistic' we mean a flow with 7 = | and 7 = f respectively, according to the 
terminology used in FKR. 




(7) 
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Fig. 1: Parameter space division for multi-transonic black hole a<;cretion. [V3] available for such accretion decreases with increase of 7, 
see text for details. 

Consider the transonic aeeretion passing through the outer sonic point rout- We take any [Ps] G C [V3] and integrate the flow for 
such [P3] upto a distance S from the event horizon and calculate the values of various V^^^ , like {M{r),T{r), p{r),p{r)}'^^^ 
etc. at rs- We developed an efficient numerical code, which will automatically pick up all [Vs] € C [Pa] out of the entire 
set [Vs] £ A[V-a], and for each [Vs] € C ['Pa], it will calculate any V^"^^ at rg for any small value of S. For S = 0.01 and 
33,1.43,1.53' Figure 3 shows the QTV of flow Mach number M(r), flow temperature T (in degree Kelvin), rest mass flow 
density p (in gm cm~^) and flow pressure p (in dyne cm~^) for shock-free accretion passing through rout- {£, A} are plotted in 
geometric unit. Any other V^'^^ for any [V3] can be obtained for S « 0.01 as well with higher computational cost. Similar 

calculations can be done for mono-transonic flows as well. For flows with shocks, general proflle of V^"^^ as a function of [Vs] 
remains fairly unaltered with the following difference of the numerical values of V^^^ for shocked and shock-free flows: 

Shock ^ No- Shocks \^ 'P ) Shock 

frpQTV QTV „QTVl 

I 'P 'P S No-Shock 
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Fig. 2: Solution topology for multi-transonic BH accretion. Sonic parameters (calculated) and the [Pa] (used) are shown at the top panel 
of the figure. See text for detail. 



We observe that the QTV of u at 0.01 Vg may become as much high as 99 percent of the velocity of Ught, leading to very 
large values of flow Mach number near the horizon. The disc thickness in the close proximity of the event horizon comes out 
to be of the order of few Kilometers. The figure is drawn for one Eddington rate accretion on a lOM© BH. Our generalized 
calculation allows us to obtain all such V^"^^s for any accretion rate onto black holes of any mass. The following trend is 

observed for the variation of V^^^ with [Pa]: 

M«^^ oc ^,u'^^'' oc ^ _^^r^QTV ^^^^^ 



A 

where oc indicates the proportionality. 

Strongly rotating purely non-relativistic flow with high energy content produces the maximum disc temperature as well as 
flow pressure at the event horizon, whereas the ultra-relativistic strongly rotating flow with lower energy content will produce 
high density ultra-fast accretion. Disc thickness is higher for purely non relativistic weakly rotating flows with large value of 
rest mass energy at inflnity, indicating the possibility of quasi-spherical inner disc structure for multiple stellar wind-driven 
accretion onto black holes. For shocked accretion, the tendency of forming the quasi-spherical inner disc region is much more 
prominent because the entropy generated at shock will puff up the post shock disc structure to provide a funnel like surface 
which supports the natural coUimation of accretion powered jets (Das & Chakrabarti 1999). 
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4 DISCUSSIONS 

4.1 Predictions for viscous flow 

In this work, viscous transport of the angular momentum is not explicitly taken into account. Even thirty years after the 
discovery of standard accretion disc theory (Shakura & Sunyaev 1973), exact modeling of viscous multi-transonic BH accretion, 
including proper heating and cooling mechanisms is still quite an arduous task, even for post-Newtonian flow, let alone for 
general relativistic accretion. Nevertheless, extremely large radial velocity close to the BH implies Tin/ << Tyisc {"Tinf and 
Tvisc are the in fall and the viscous time scales respectively), hence our assumption of inviscid flow is not at all unjustiflcd, at 
least upto a few to tens of r-g or so. Far away from the BH, this may not be a very good assumption. However, one understands 
that one of the most significant effects of the introduction of viscosity would be the reduction of radial angular momentum. 
We found that the location of the sonic points anti-correlates with A (weakly rotating flow makes the dynamical velocity 
gradient steeper), which indicates that for viscous flow the sonic points will be pushed further out and the flow would become 
supersonic at a larger distance for the same set of other initial boundary conditions. The terminal values of the Mach number 
and disc height anti-correlate, while the density and pressure correlates with A. Our calculations, if applied for viscous flow, 
will thus produce a more supersonic and quasi-spherical structure near the event horizon, while the terminal flow pressure 
and density would be reduced. 



4.2 Spectral Signature of Black Hole Spin 

One of the most tantalizing issues in black hole astrophysics is to study whether the spin of the black hole (the Kerr parameter 
a) could be determined using any observational means, the spectral signature of the black hole, for example. Our investigation 
of general relativistic accretion flows in Kerr metric may provide an important step towards a better understanding of this 
problem. For accretion onto the Schwaxzschild black holes, we have been able to track the infalling matter upto the point 
extremely close to the event horizon. We can successfully perform the same procedure (with increasing degree of mathematical 
complexity) for accretion onto Kerr black hole as well and can calculate various dynamical and thermodynamic flow variables 
upto the point extremely close to the event horizon (as close as it could be) as a function of Kerr parameter a (T. K. Das, 
Paul J. Wiita & Paramita Barai, in preparation). Hence the terminal behaviour of infalling matter (e.g., flow temperature, 
pressure, density etc. which are responsible to characterize the observed spectra) could be studied as a function of the black 
hole spin and the spectral signature of the black hole rotation may, at least qualitatively, be determined. This, will, as we 
believe, be an important step forward towards a better understanding of black hole astrophysics. However, we would like to 
mention here that such spectral signatures would be very difflcult to detect observationally. Most of the radiations produced 
such a close vicinity of the event horizon will be directly swallowed by the black hole itself, rest of the radiative signature will 
eventually suffer an enormous amount of gravitational redshift, and hence will remain almost inaccessible for analysis using 
any present day observational techniques. Also, the temperature profile obtained from our work will mainly correspond to the 
extreme ultra-violet or X-rays, while most of the observed flux for black hole candidates come from flares and jets instead of 
the inner region of the disc, and hence the theoretically calculated strong gravity effects would be diluted. 



4.3 Hawking radiation from acoustic black holes 

The Pioneering attempt made by Unruh (1981) on mapping certain aspects of BH physics onto the propagation of sound wave 
in supersonic convergent flows leads to a number of important work (see Visser 1998, V98 hereafter, and references therein) to 
identify the propagation of acoustic disturbance in a polytropic inviscid fluid flow with the d'Alembertian equation of motion 
of a minimally coupled massless scalar field propagating in a (3-1-1) Lorentzian manifold. Sound propagation is described 
by an acoustic metric algebraically dependent on the flow density and Mach number; which is conformaly related to the 
Painleve -Gullstrand-Lemaitre representation of Schwaxzschild geometry, with a constant conformal factor close to the black 
hole event horizon. Such an equivalence is potentially useful to tackle the tantalizing issue of Hawking radiation in terms of 
physical quantities governing the transonic fluid flow. In this paper we construct multiple sonic surfaces (a collection of r^n and 
Tout forming a 3D hypersurface on \P-i\ space) which is essential for a complete investigation of the formation of outer-trapped 
surfaces in acoustic ergo-region and of acoustic horizons. We investigate the multi-transonic, barotropic inviscid fluid flow 
(accretion) using a complete general relativistic framework and not using any Newtonian assumption. Hence, we believe that 
our work presented in this letter, may have a two fold importance. Firstly, a general relativistic calculation of the values of 
the fundamental dynamical and thermodynamic flow variables close to the black hole event horizon will help in constructing 
a diagnostic high energy spectra, useful for the identification of galactic and extragalactic black hole candidates. Secondly, a 
full general relativistic treatment of representative transonic fiuid flow will take an welcome step-forward towards a concrete 
formulation of the 'Acoustic General Relativity', eis the term used by V98. We have also established a very important fact 
that general relativistic spherical accretion onto astrophysical black holes may generate acoustic white holes, and the ratio 
of the Hawking temperature to the analogous sonic black hole temperature is independent of the black hole mass, which is a 
nice way out from the great debate whether astrophysical/ primordial black holes can really have such a lower mass limit so 
that the black hole micro-physics can really be dealt with. Details of such calculations are beyond the scope of this paper and 
have been presented elsewhere. 
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